Coexistence of Spin Density Wave and Triplet Superconductivity 
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We discuss the possibility of coexistence of spin density wave (antiferromagnetism) and triplet 
superconductivity as a particular example of a broad class of systems where the interplay of mag- 
netism and superconductivity is important. We focus on the case of quasi-one-dimensional metals, 
where it is known experimentally that antiferromagnetism is in close proximity to triplet supercon- 
ductivity in the temperature versus pressure phase diagram. Over a narrow range of pressures, we 
propose an intermediate non-uniform phase consisting of alternating antiferromagnetic and triplet 
superconducting stripes. Within the non-uniform phase there are also changes between two and 
three dimensional behavior. 

PACS numbers: 74.70.Kn, 74.25.DW 
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The competition or coexistence of magnetic order and 
superconductivity is a very important problem in con- 
densed matter physics. There is a broad class of systems 
that present magnetic order and superconductivity in 
close vicinity. One of the most important systems are the 
Copper Oxides, where singlet superconductivity is found 
next to antiferromagnetism In addition, striped 

phases, where coexistence of antiferromagnetic order and 
singlet d-wave superconductivity, were observed in Cop- 
per Oxides 0- Another system where magnetism and 
superconductivity are intertwined is Strontium Ruthen- 
ate, where the proximity to ferromagnetism has been ar- 
gued as being important to the existence of triplet su- 
perconductivity in these materials Furthermore the 
newly discovered ferromagnetic superconductors ZrZri2 
and UGe2 have stimulated a debate on the coexistence 
of ferromagnetism and triplet or singlet superconductiv- 
ity 0, . However, unlike any of these previous exam- 
ples, we will discuss in this manuscript a system which 
may exhibit coexistence of antiferromagnetism (AF) and 
triplet superconductivity (TS). 

New experiments on quasi-one-dimensional supercon- 
ductors in high magnetic fields have shown that TS is 
strongly affected by the proximity to an AF phase char- 
acterized by insulating spin density wave (SDW) or- 
der 0. From now on we will use interchangeably SDW 
and AF. Motivated by these experiments and the known 
phase diagram of quasi-one-dimensional (TMTSF^PFg 
under pressure we propose a new phase for quasi-one- 
dimensional systems where AF (SDW) and TS coexist. 
The coexistence of these phases implies that the new 
state is non-uniform, with alternating stripes of SDW 
and TS, due to the appearance of a negative interface 
energy between SDW and TS regions. As indicated in 
the schematic phase diagram (Fig.l), the inhomogeneous 
intermediate phase is expected to exist over a narrow 
range of pressures AP = P2(T)—Pi(T) around P Cl where 
AP <c P c . 

Effective Free Energy: The possibility of coexistence 
of SDW and TS in quasi-one-dimensional conductors 
transcends microscopic descriptions based on standard 



g-ology, where SDW and TS phase boundaries neigh- 
bor each other but do not coexist Inspired by ex- 
periments HII(3> we model (TMTSF) 2 PF 6 as a highly 
anisotropic orthorombic crystal, and we take the primary 
directions of the SDW vector order parameter to be the 
b-axis (y-direction), and the primary direction of the TS 
vector order parameter to be the c-axis (z-direction). 
Furthermore, we consider the spatial variation of the 
SDW or TS order parameter to be along the a-axis (x- 
direction), as a reflection of the quasi-one-dimensionality 
of the system. This simplifies the choice of the order pa- 
rameters to be S(r) — > Sb(x), and D(r) — > D c (x), and 
reduces the associated effective field theory to one spatial 
dimension. Thus, the generalized Ginzburg-Landau Free 
energy in real space can be written as 



J 7 tot = -F AF + J~TS + Fc, 



(1) 



where Faf, 3~ts, and Tc are the AF, TS and coupling 
contributions discussed below. The AF contribution is 



•F af = 



dx [Uaf(x) + V A f{x)} , 



(2) 



Lab 



where Uaf(x) = c*af 

1af\ Sbix)] 4 represents a typical GL Free energy den- 
sity, and Vaf = S A F\d x S b (x)\ + 9 AF \S b (x)\ \d x S b (x)\ 
represents the extra terms in the expansion, which are 
relevant close to Pi(T) (Fig. 1). The TS contribution is 



dx [Uts(x) + V TS (x)} 



(3) 



Lts 



where Uts{x) = olts \D c {x)Y + (3 T s\d x D c (x)\ z + 
1ts\ D c (x) 1 4 represents a typical GL Free energy den- 
sity, and V TS = S TS \d x D c (x)f + 6 TS \D c {x)\ 2 \d x D c {x)\ 2 
represents the extra terms in the expansion, which are 
relevant close to PiiT) (Fig. 1). To describe the coexis- 
tence region the two order parameters must couple. To 
conform with independent Parity invariance, 



?C=J2 I dx\' bc \S b (x)\ 2 \D c (x)\ 2 



(4) 
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FIG. 1: a) Phase diagram of (TMTSF) 2 PF 6 in a log-linear 
scale (from [3), showing schematically the proposed SDW- 
TS coexistence region (inset), b) Schematic drawing of the 
proposed stripe pattern for the SDW-TS coexistence region. 

where the sum is over all interfaces between AF and TS, 
the coupling constant X' bc is pressure and temperature 
dependent, and £ p is the proximity length over which AF 
and TS order parameters coexist locally. This length can 
be written as £ p — £ Pt AF + @p,ts > where £ Pi af is the AF 
proximity length into the TS region, and £ p ,ts is the TS 
proximity length into the AF region (See Fig. 2). If 
X' bc > it is more favorable for the AF and TS phases 
to phase-separate, however if X' bc < an inhomogeneous 
phase with a large number of interfaces is favored, and 
the coexistence of AF and TS is possible. 

In non-triplet systems local AF order and local sin- 
glet superconductivity (SS) can in principle coexist since 
AF order favors singlet correlations, and the proximity 
lengths on SS/AF systems can be small or large de- 
pending on the SS and AF materials 0, 0. How- 
ever, we are interested in TS and not in SS. In this 
case, it is well known that AF order is pair breaking 
to triplet electron pairs an d it is expected that 

£p,AF and £ p ,ts are small in comparison to the lengths 
of £ts and £af of the TS and AF stripes, respectively. 
Only when the pressure P is close to the phase bound- 
aries Pi or P 2 (shown in Fig. 1) where £ts and £af 
approach zero respectively, the proximity lengths £ p .ts 
and £ p af can be comparable to £af and £ts- Using 
WKB [TlJ and the deGennes extrapolation methods 



the upper bound for i p is £ P (P) < 0.1[£ AF {P) + £ts(P)}, 
and £ p ,af{P) < O.U TS (P 2 ) and £ p ,ts{P) < 0.1£ AF (Pi)- 
Thus, extrapolating AF and TS order parameters in the 
proximity region by linear functions a la deGennes [lij 
leads to 

Tc = hc\d x S b (x)\ 2 \d x D c (x)\ 2 , (5) 

inter 

where Xb c — ^'bcIo P dxx 2 (£ p — x) 2 is the new coupling 
constant. This coupling describes well the inhomoge- 
neous phase where proximity effects between AF and TS 
stripes are weak, i.e., away but not too far from phase 
boundaries P\ and P 2 . 

Proximity effects will be important close to either 
phase boundary (Pi,P 2 ) as they may lead to further cou- 
pling between like-stripes. The inter-TS-stripe coupling 
(Josephson-type) 

Fi,ts = ^ / dxri T s\D c . n+1 (x) - D Ctn (x)\ 2 (6) 

n J overlap 

is significant for P ~ P 2 , where the system changes from 
2D to 3D TS. The inter- AF-stripe coupling 

Fi,af = X/ / dxriAF\Sb, n +i(x) - S b ^ n (x)\ 2 (7) 

n J overlap 

is significant for P « Pi, where the system changes from 
3D to 2D AF. The domain of integration for both cases 
above is the overlap region between two consecutive AF 
or TS stripes respectively. 

Saddle Point Equations: To obtain the saddle point 
equations, we consider J-q, J'i^af and Fi,ts perturba- 
tively and minimize J- to t with respect to Sb(x) and D*(x). 
Variations of Ttot with respect to Sb(x) lead to the dif- 
ferential equation 

[2a AF + A lAF Sl{x) ~ l3 A Fd 2 ]S b (x) + M AF S b {x) = 0, 

with M A FSb(x) = -SafOx [(d x S b (x)] 3 + 
20AFSb(x)\d x Sb{x)\ 2 . Variations of JT tot with re- 
spect to D*(x) lead to a similar equation. In the case 
where X' bc {P,T) < 0, the formation of an inhomogeneous 
phase of alternating AF and TS stripes is preferred, and 
two additional transition lines (Pi, P 2 ) emanate from 
(P C ,T C ). The presence of the inter-TS and inter- AF 
stripe Free energies indicates that D n+ i{x) — D 71 (x+xq), 
and S n +i{x) = S n (x + xq) (in phase solutions), since 
rjAF and rjTs are both positive. For such inhomogeneous 
phase, the boundary conditions in the presence of AF-TS 
interfaces can be chosen as in deGennes method 0] by 
requiring that S b (x)\ mter + = and D c (x)\ mter - = 0, 
where inter + and inter~ denote the two boundaries 
limiting the region of locally coexisting Sb(x) (AF) and 
D c (x) (TS) (See Fig. 2). 

Variational Free Energy: We analyse J-tot variation- 
ally. We consider first the AF case and search for pe- 
riodic solutions with period £af, with Sb(x)\ inter + = 
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FIG. 2: Alternating AF and TS stripes with local coexistence 
region £ p , which is defined by the region between interface + 
and interface ~ where AF and TS order parameters vanish. 
The length £ p = £ p> af + £ p ,ts, where £ p ,af and £ p ,ts are AF 
and TS proximity lengths. 



at the AF-TS interfaces. For a given volume of the AF 
region, controlled by Laf, the Free energy associated 
with the AF phase becomes the sum of Naf identical 
terms, where Naf — Laf/^af gives the number of AF 
stripes. Generally, each term in Taf corresponds to an 
insulating AF stripe characterized by the order parame- 
ter Sb(x) = J2n A n sin(Q„2;), where Q n = 2-Kn/i A F- But 
here we take for simplicity the variational class where 
Sb{x) = A\ sin(Qix). To simplify notation, we will just 
use A\ — > A and Qi — > Q. In this case 



Taf — Lai 



C 2 (Q)A 2 + C 4 (Q)A 4 



(8) 



where Cfe(Q) = (a AF + PafQ 2 )/2 and C 4 (Q) = (3j AF + 
@afQ 2 + 3SafQ 4 )/&- The same type of analysis applies 
to Tts- In the absence of magnetic field, we assume pe- 
riodic solutions of the form D c (x) = Y] B n sm(K n x). 
As in the AF case we confine ourselves to a single com- 
ponent variational form D c (x) — Bism(Kix), and use 
the simplifying notation B\ — > B and K\ — > K. Here B 
can be complex, but independent of position x. All the 
analysis discussed for the AF (SDW) case applies with 
the following change of notations: Laf — * Lts, A — » B, 
Q — » K, aAF — > o>ts, Paf — > (3ts, etc. Which in the TS 
case leads to: 



T TS = L TS D 2 (K)\B\ 2 + D 4 (K)\B\ 



(9) 



where D 2 (K) = (a TS + P T sK 2 )/2 and D 4 (K) = (3~/ T s + 
OtsK 2 + 36tsK 4 )/8. And the coupling Free energy is 

T C =N mt A(Q,K,£ p )A 2 \B\ 2 , (10) 

where Ni nt = 2N is the total number of interfaces, fi n t = 
A(Q, K, £ P )A 2 \B\ 2 is the Free energy of one interface with 
A(Q,K,£ P ) = X' bc J e dx\sin(Qx) sin[K(xQ +a;)]| 2 , where 

£o = £ts ~ V 

Variational Solution: Variations of Ttot with respect 
to 4>af — A, tpTs = \B\, and qAF = Q or qxs — K lead 
to the non-trivial solutions 

, 2 4ft0; - 2Aa i 8 i 



(11) 



a, 



ft0, - 6aA 
The width of each stripe then is given by 



L = 2tt, 



! Pi' 



6a, Si 



a; 



(12) 



(13) 



where i = AF, TS. In the case of X' bc < 0, the transition 
line Pi (T) corresponds to the disappearance of the pure 
AF (SDW) phase, and the transition line P 2 (T) corre- 
sponds to the appearance of the pure TS phase. This im- 
plies that at Pi(T) the TS stripe width is £ts = 0, while 
at P 2 (T) the AF stripe width is £af — 0. Furthermore, 
for Pi(T) < P < P 2 (T), £ TS (£af) increases (decreases) 
with increasing pressure. In order to meet these and the 
saddle point requirements, the parameters appearing in 
Ttot must behave as follows. We define the reduced pres- 
sure changes AP m = [P — P m (T)]/P c , where m = 1,2 
and the density of states N(Ep) at the Fermi energy Ep 
to analyse the AF and TS parameters. For P < P 2 (T), 
the AF parameters have the form jaf — 1iN(€f)T 2 > 



6iN(e F )T t 



2 > 0; 

21 



VAF 



ViN(e F )T 2 > 



0; 5 AF 

0; a AF = aiN(e F )T 2 \AP 2 \ e <*AF , with ai < n. 
Paf = ft^(e F )T 2 |AP 2 | e ^, with ft < 0; 9 AF 



e l N{e F )T 2 \AP 2 \ E <>AF , with 9i < 0; and 36j AF 5 AF - 
9\ F > 0. For P > Pi(T), the TS parameters have the 
form 7 Ts = l2N{e F )T 2 > 0; S TS = 6 2 N{e F )T 2 > 0; 
Vts = V2N(e F )T 2 > 0; a TS - a 2 iV(e F )T c 2 |AP 1 | e «T S , 
with a 2 < 0; (3 TS = f3 2 N{e F )T 2 \APi\ e ^T S , with 



ft < 0; 9 TS = 6» 2 A^(e F )r c 2 |AP 1 | £( 'T S , with 6> 2 < 0; 
and 36jts$ts ~ @rs ^ ®- Consider now, the in- 
terface terms in the region P\{T) < P < P 2 (T), 
which has the form X' bc = X N(e F )T 2 sgn[(P - P X )(P - 
P 2 )]|APi| e *Af |AP 2 | £a ts, with A > 0. This form is re- 
quired to make the interface energy negative between 
P(T) and P 2 (T). 

Phase Boundaries: Next, we focus only on the 
analysis of £af and Ttot in the vicinity of P 2 {T). We 
note in passing that the analysis of £ts and Ttot in 
the vicinity of Pi(T) is entirely analogous. Under 
these considerations, as P — > P 2 (T) the size of the AF 
stripes is given by £af ~ Wi\A.P 2 \^ Eo: af- e Paf)/ 2 , where 
W x = 27r(cv 1( 5i/ft7i) 1/2 , when e 0AP + e 0AF > e aAF , and 
£a A F + £«af > £ Paf- The requirement that — > as 
P — > P 2 (T) forces £ QA f > £ ^f ■ Since the number of AF 
and TS stripes are the same Naf — Nts = N, where 
N = L/(£ AF (P) + £ts(P) - i P (P)), the number of in- 
terfaces is Nmt = 2N. Therefore, the four contributions 
to Ttot in the coexistence region are Taf = Laf/af, 
Tts = Lrsfrs, Fc = 2Nf int , Tprs = 2Nfi t TS with 
L AF {p) = N£ AF {P), L ts = N£ TS {P). As P -> P 2 , 
the AF stripe length £af{P) —> 0, while the TS stripe 
length £ T s{P) -> ^ts(P2(T)). Let us analyse the 
behavior of Ttot near P 2 (T) term by term. The AF part 
Taf is the product of L AF (P) w L^ F (P)/^ TS (P 2 ), 



where £af(P) 



Const.\AP 2 \ ( - e "AF-^AF)/ 2 l and 
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FIG. 3: a) Normalized amplitudes of the order parameters 
(A = A/\A(Pi)\, B = B/\B{P 2 )\) and stripe lengths Iaf 
and Its normalized by £af(Pi) and Its(P2), for exponents 
e ai — 5.0, Ep i — 1.0, Eo i — 5.0, where i = AF, TS and 
dimensionless parameters: Si = /3i = 1.0, 71 = 0.07, 5i = 
0.007, 6>i = -0.003; and & 2 = f3 2 = 1.0, 72 = 0.06, S 2 = 
0.007, 9 2 — —0.002. b) Free energies for the coexistence and 
the pure AF and TS phases for the same parameters of a) 
and £ p = 0.05(^af + Its), e Xi = 1-0, A = -0.003, and 
f)i — fj 2 — 2.9 x 10~ 6 . Solid dots — + Ttot; stars — > Lfrs; 
triangles — > LJaf', squares — > Tc- P[ indicates change from 
3D to 2D AF, and P 2 denotes change from 2D to 3D TS. 

f AF = -Const.\AP 2 \ 2£ ' 3 AF , which leads to 
Taf = -Const.\AP 2 \ {ea AF+^fiAF')/ 2 . The TS part 
T TS is the product of L TS (P) « L[l - Iaf{P) / Its{Pz)] 
and f TS (P) « /ts(P 2 ) < 0, thus .F TS » 
Fts(P2) + LConst.\AP 2 \ {e °'AF-^AF)/ 2 i where 
•7"ts (-P2) = LJts{P2) is the Free energy of the 
pure triplet phase. Thus, for P < P 2 , Tts is increased 
with respect to the pure phase Tts(P2)- Furthermore, 
Tc = —Const.\P 2 — P\ £x af+ £ Paf is negative, and 
Fi,ts is positive but proportional to a higher power 
of AP 2 . Taf can be neglected in the vicinity of P 2 
because it depends on a higher power of IAP2I than 
that of Tts or Tc, when e aAF > ep AF . For T to t to 
be lower than that of the pure TS phase Tts(P) it is 
necessary that the negative interface Free energy Tc 
dominates. This imposes the following requirement 
£ aAF > 3 £ /3af + 2e\ AF . If ex AF + ep AF > 1 the phase 
transition at P 2 is continuous. If e\ AF + £p AF < 1 
the phase transition at P 2 is discontinuous. Similar 



analysis close to P\ leads to a continuous transition 
when £\ TS + £[3 TS > 1 and to a discontinuous transition 
when e\ TS + sp TS < 1. Thus, the point (P C ,T C ) can 
be bicritical, tricritical or tetracritical. In Fig. 3, we 
show the behavior of the various contributions to Ttot 
for the case where the transitions are continuous at 
Pi and P 2 , and (P C ,T C ) is tetracritical. Dimensionless 
"volume" parameters are defined as en — fli = p 1 / 2 , 
7i = Ji^i 1/2 , ^ = £ 4 of /2 , 0i = 6iv] 12 , "surface" 
ones are defined as Ao = Aqctj 1 ^ 4 ct 2 1 ^ 4 (^ p (Pi)ct^ 2 ) 5 , 

% = ^^^^^(^(POcry 2 ) 3 , where pi = Oifa, a t = a i /p i , 
and i = l,2. 

Final Comments: A more realistic description of the 
system should include variations of SDW and TS order 
parameters along the y and z (transverse) directions, and 
be highly anisotropic but truly three dimensional. This 
is important for a renormalization group (RG) analysis 
and the determination of critical exponents in the case 
of continuous transitions. 

Summary: We have proposed the possibility of coexis- 
tence of antiferromagnetism and triplet superconductiv- 
ity in the phase diagram of (TMTSF) 2 PF 6 . This inter- 
mediate phase is proposed to be inhomogeneous and to 
consist of alternating insulating AF and TS stripes 01 
Two additional transition lines appear in a narrow range 
of pressures around P c separating the coexistence region 
from the pure AF and pure TS phases. We estimate the 
maximum pressure range to be AP/P C w 10% at T = 0. 
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